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Abstract 

Background: Antibiotic resistance in bacterial infections is a growing threat to public health. Recent evidence 
shows that when exposed to stressful conditions, some bacteria perform higher rates of horizontal gene transfer 
and mutation, and thus acquire antibiotic resistance more rapidly. 

Methods: We incorporate this new notion into a mathematical model for the emergence of antibiotic multi- 
resistance in a hospital setting. 

Results: We show that when stress has a considerable effect on genetic variation, the emergence of antibiotic 
resistance is dramatically affected. A strategy in which patients receive a combination of antibiotics (combining) is 
expected to facilitate the emergence of multi-resistant bacteria when genetic variation is stress-induced. The 
preference between a strategy in which one of two effective drugs is assigned randomly to each patient (mixing), 
and a strategy where only one drug is administered for a specific period of time (cycling) is determined by the 
resistance acquisition mechanisms. We discuss several features of the mechanisms by which stress affects variation 
and predict the conditions for success of different antibiotic treatment strategies. 

Conclusions: These findings should encourage research on the mechanisms of stress-induced genetic variation 
and establish the importance of incorporating data about these mechanisms when considering antibiotic 
treatment strategies. 

Keywords: stress induced mutagenesis, HGT, antibiotic resistance, evolution, mathematical model 



Background 

Bacterial resistance to antibiotics has accompanied the 
introduction of new antibiotics since shortly after peni- 
cillin was first introduced [1] and is currently considered 
a major health issue [2,3]- The implications of infection 
with antibiotic-resistant bacteria include increased mor- 
tality rates, prolonged hospitalization and higher cost of 
treatment [1,4,5]. A particularly dangerous prospect of 
the continued evolution of drug resistance in bacteria is 
the creation of new, multidrug resistant bacteria. Such 
bacterial strains are already present in several species of 
bacteria [3,6] and treating them is more difficult and 
often accompanied by a period of ineffective treatment, 
resulting in increased patient mortality [6]. Moreover, 
the rate of new drug development is declining, leaving 
few treatment alternatives for treating the increasing 



* Correspondence: lilach.hadany@gmail.com 

Department of Molecular Biology and Ecology of Plants, Tel Aviv University, 
Ramat Aviv, Tel Aviv 69978, Israel 



burden of multi-resistant bacteria [7,8]. Since resistance 
is especially prevalent in hospitals [9], various treatment 
strategies have been suggested to facilitate better 
responses to resistant infections and minimize the emer- 
gence of new multi-resistant bacteria. 

Three prominent strategies of antibiotic treatment are 
cycling, mixing and combining. Under a cycling regime, 
all the patients are treated with the same antibiotic drug 
at a given time, and the drug used is periodically 
switched. The rationale behind cycling is that each time 
an alteration of drugs is administered, the pathogens 
resistant to the previously used drug are attacked and 
are hopefully susceptible to the new drug [10]. In the 
mixing strategy, each patient receives a randomly 
selected drug. This strategy can be viewed as the default 
antibiotic usage within a hospital unit, when there is no 
preference for any particular antibiotic. In such a case, if 
two relevant antibiotics exist, approximately half the 
patients receive each drug at any given time. Mixing has 
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the advantage of creating a heterogeneous stress envir- 
onment for the bacterial population [10]. At each trans- 
mission, a bacterium has a probability of half 
encountering a drug to which it has not recently been 
exposed, and hence to which it is unlikely to be resis- 
tant. Combining is the administration of several drugs 
to each patient. By applying several antibiotics at once, 
combining is designed to diminish the chance of evol- 
ving resistance by eradicating any bacteria resistant to 
just one type of antibiotics. As a result, more antibiotics 
are used in combining than in mixing or cycling. This 
could lead to higher antibiotic-related toxicity and 
increased treatment costs [11]. 

Attempts to compare the different treatment strategies 
and assess their relative efficiency have been made both 
in empirical studies [12-17] and using theoretical analysis 
[10,18-20]. All in all, results obtained using both 
approaches have been inconclusive. It seems clear that 
current models do not capture all the aspects of the phe- 
nomenon, thus failing to properly distinguish between 
scenarios favoring different treatment strategies. We sug- 
gest that part of this shortcoming may result from the 
simplifying assumption that resistance is acquired at a 
constant rate, ignoring recent evidence to the effect of 
environmental stress on the mechanisms of resistance 
acquisition. The frequency of horizontal gene transfer 
(HGT) and mutation was shown to increase when bac- 
teria are under various stressors, including nutritional 
deprivation, DNA damage, temperature shift, oxidative 
stress and exposure to antibiotics [21-24]. An environ- 
mental stressor especially relevant to our subject of 
inquiry is antibiotics. In Streptococcus uberis, acquisition 
of rifampin resistance through mutation was shown to 
increase more than 1,000-fold under ciprofloxacin [25]. 
Interestingly, the rifampin resistant mutants showed no 
resistance for ciprofloxacin, which indicates that this was 
not merely the result of selection. It was also shown that 
Pseudomonas aeruginosa increases its mutation rate by 
up to 10 in the presence of tetracycline antibiotics, and 
consequently obtains resistance to antibiotics [26]. 
Stress-induced mutation (SIM) might therefore have a 
substantial influence on the dynamics of antibiotic resis- 
tance acquisition [27-29]. In the context of HGT, Strepto- 
coccus pneumoniae was shown to increase the rate of 
chromosomal DNA uptake by transformation of a mar- 
ker conferring resistance to streptomycin, when treated 
with either streptomycin or norfloxacin [30] . It was also 
shown that ciprofloxacin induces the transfer of the SXT 
integrating conjugative element, which is known to 
encode for antibiotic resistance genes, in SXT-containing 
Escherichia coli and in Vibrio cholerae up to 300-fold 
[31]. Phages were also observed to increase horizontal 
transfer of genetic material as a reaction to their host's 
antibiotic-induced SOS response[32], a process which 



might lead to an increased rate of antibiotic resistance 
acquisition [33]. Theoretical work also supports stress- 
induced genetic variation as a successful evolutionary 
strategy [34-36] so this phenomenon might be even more 
widespread then we currently know. Our goal is to 
explore the impact of genetic variation induced by anti- 
biotic stress on the spread of antibiotic multi-resistance 
in a hospital unit. We use a classical modeling approach 
(first described in [37]), modified to describe SIM and 
stress-induced HGT. Our model is used to evaluate the 
efficacy of each treatment strategy under different 
assumptions regarding the effect of stress on genetic var- 
iation. We find that stress-induced variation can indeed 
alter the preferred treatment strategy. 

Methods 

Our mathematical model describes the dynamics of bac- 
terial infections in a hospital unit. The bacterial pathogens 
in question are assumed to accompany other ailments and 
not be the main reason for hospitalization. We consider 
two different antibiotic drugs, denoted antibiotic 1 and 
antibiotic 2. The frequencies of patients infected with bac- 
teria resistant to antibiotics 1 and 2 are 7? x and R2 , respec- 
tively, and the frequency of patients infected by 
susceptible bacteria is S. The frequency of uninfected 
patients is X . Clearance due to antibiotic usage occurs at 
rate x , and X; determines the fraction of patients receiving 
antibiotic i. Resistance is assumed to be complete, so that 
a patient infected by a bacterial strain resistant to drug 1 
will not be affected at all by treatment with that drug. 
Conversely, if treated with drug 2, the patient becomes 
uninfected ( X ) at rate x . y is the rate of spontaneous 
clearance due to the response of the patient's immune sys- 
tem, [} is the rate of bacterial transmissions resulting in 
infection (for simplicity, superinfection is neglected), and 
m is the rate of patient turnover (so that the mean patient 

hospitalization time is — days). Since we assume the bac- 
m 

terial infection is not the main reason for hospitalization, 
patients leave the hospital or die at a rate proportional to 
their frequency. The proportion of infected patients enter- 
ing the hospital is determined by X^, Xr 2 , Xs for patients 
carrying bacteria resistant to antibiotic 1, to antibiotic 2, 
or to none, respectively. Uninfected patients enter the hos- 
pital at rate m (l — (Xr 1 + Xr 2 + Xs)), so hospital occu- 
pancy is kept constant. All of the parameters representing 
rates are given in units of day" . We assume there are no 
double-resistant bacteria in the hospital initially, and that 
their frequency in the general population is negligible. 
This scenario may reflect situations where newly devel- 
oped antimicrobial agents have been recently introduced, 
or were kept as the last resort, so that double resistance is 
still scant. 
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These parameters are incorporated in the following set 
of ordinary differential equations: 

— = (1 - X s - X Rl - Xr 2 - X)m + (t + y )S + (t/2 + Y )R\ 
at 

+{txi + y)R2-PX{S + R 1+ R2) 

■^ = (ks-S)m-{T + r)S + fSXS (E1) 

dR 
dt 

= - Ri)m - (rxi + Y)Ri + PXR 2 

at 

The equations describe the rate of change of patient 
frequencies within a hospital. The dynamics are illu- 
strated in Figure 1. 



{k s - S)m - {x + y)S + PXS 
- = (A.r, - Ri)m - (t/ 2 + Y )Ri + PXRr 



Equations El were solved analytically [See Additional 
file 1, section 3], and by numerical integration using 
Matlab® R2009a. 

Moving average calculation: We used the numerical 
solutions of equations El and the analytical computations 
of double resistance emergence [See Additional file 1, 
sections 1 and 2] to calculate values of a moving average. 
First, the different sets of parameters are ordered accord- 
ing to the parameter plotted on the x axis. Each point in 
the plot presents the average double resistance emer- 
gence over 201 equally weighted parameter sets: the one 
corresponding to the point itself, the 100 nearest para- 
meter sets with lower values of r(— — ), and the 100 
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Figure 1 Illustration of antibiotic resistance dynamics in a hospital setting. The solid lines represent infection, recovery and patient 
turnover. Dashed lines represent the effects of HGT and mutation. Stress-induced genetic variation would result in an increased weight of the 
dashed lines under antibiotic stress. X, 5, fl q , R 2 are the frequencies of uninfected patients, patients infected with susceptible bacteria, and 
patients with bacteria resistant to antibiotics 1 and 2, respectively. They enter the hospital with rates 

(1 — A. s — Xr, — Xn 2 )m, X s m, Xr 1 m, Xn 2 m, and leave with a rate proportional to their frequency, where the patient turnover rate, m , is 
the proportion constant. Infected patients turn uninfected either through spontaneous recovery due to the immune system (at rate y), or due to 
antibiotic treatment (at rate r). and %i determine which amount of antibiotics 1 and 2 are used, respectively. Uninfected patients become 
infected at rate /3 multiplied by the frequencies of cleared and infected patients. /? 12 denotes the fraction of patients infected with double 
resistant bacteria, assumed to be zero at the beginning. HGT, horizontal gene transfer. 
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nearest parameter sets with higher values of r( — -) . The 

same 201 data sets are used to calculate double resistance 

emergence for a given value of x ( — — ) under each of the 

three strategies, resulting in a correlation between the 
three values, and a similarity in the shape of the three 
moving average curves (since the three means are taken 
over the same parameter sets). 

Results 

Stress-induced mutation 

Considering stress-induced mutation, we define fi s and 
fi r as the rates of mutations conferring antibiotic resis- 
tance when bacteria are under antibiotic stress and 
when they are free of antibiotic stress, respectively. We 
assume that mutation is increased by stress, that is, fA s 
>ft r . Our work concentrates on the scenario where sen- 
sitive bacteria (corresponding to host frequency S ) have 
a negligible contribution to the generation of single 
resistant bacteria. In terms of our dynamics, this trans- 
lates to the following inequality: 

U, S S « X Rl m + pXRi, i=l,2 (CI) 

Intuitively, inequality CI is satisfied due to either 
abundance of single resistance in the population outside 
the hospital, causing high entrance rates of single resis- 
tant bacteria, or abundance inside the hospital due to 
infection and selection. 

An important feature of antibiotic resistance is its per- 
sistence within a host without direct selective forces for 
long periods of time [38]. Direct selection occurs when 
a patient is treated with a certain antibiotic, say antibio- 
tic 1, and a new bacterium resistant to antibiotic 1 arises 
within the host, by mutation or HGT. Thus, due to 
strong selection for antibiotic resistance, it has a high 
probability of taking over the entire infecting population 
and turning the host to R^. When the host is not treated 
with antibiotic 1, resistance to antibiotic 1 might not 
confer any direct fitness advantage. Thus, we assume 
that the probability of a bacterium to take over the 
infection in the second scenario is a times the chance in 
the first scenario, where 0 < a < 1 represents the relative 
persistence of antibiotic resistance when there is no 
direct antibiotic usage. Under condition CI, we can con- 
sider a only when computing within-host dynamics of 
double resistant bacteria. 

The exact value of <7 is hard to measure, but is likely 
to be non-zero, as evidence suggests that patients who 
have not been treated with antibiotics for periods of up 
to three years still carry antibiotic resistant bacteria [6] . 
One reason for high persistence of resistant bacteria in 
ambulatory patients and medical staff might be residuals 



of antibiotics that are found in the environment at 
amounts sufficient to change the fitness of sensitive bac- 
teria. This might often be the case in hospitals, as it was 
shown that even very low concentrations of antibiotics 
can select for antibiotic resistant bacteria [39] and that 
even ambulatory patients who have not received antibio- 
tics for long periods of time harbor high frequencies of 
antibiotic resistant bacteria [40]. 

To compute the emergence of double resistance 
through mutation we first define a term describing the 
sum of the frequencies of patients carrying single resis- 
tant bacteria in the hospital for a certain time period: 

d 

Es,m{U) := j {R 1 {U)+R 2 (U))dt 

to 

Using E SIM and the parameters described above, we 
can now write a term describing the emergence of dou- 
ble-resistant bacteria under treatment strategy U. It will 
be denoted by C, sm {U ). Writing c, SIM explicitly for the 
three strategies results in the following expressions [See 
Additional file 1, section 2, equations 7, 8, 10]: 

ismimixing) = j{/x s + n r a)E S iM{mixing), 

%sim {combining) = h s E S im {combining) 

For cycling we have a more complex expression. We 
will divide time into segments in which only one anti- 
biotic is applied. In each of these segments only one 
strain of resistant bacteria is under antibiotic stress. 
Thus, 

£sim {cycling ) = p s A + fi r o B where A = J R t {cycling) 
when antibiotic * i is applied and B = f R t {cycling) 
when antibiotic i is applied [See Additional file 1, sec- 
tion 2, definition (9)]. This allows us to prove that as 

Ms 

— increases, emergence of double resistance in cycling 
fi r 

decreases relative to the emergence of double resistance 
in mixing and combining [See Additional file 1, section 
2.1, equations 11-12]. This is shown in Figure 2, where 
emergence of double resistance by mutation for each 

Ms 

strategy is plotted as a function of — . Dashed and solid 

Mr 

lines correspond to a = 0 and a = 1. Double resistance 
emergence is monotonically increasing in a [See Addi- 
tional file 1 section 2, equations 7, 10], thus the values 
for intermediate values of <7 are between the two curves. 
Some intuition about the relative success of cycling can 
be obtained from looking at the fraction of stressed bac- 
teria under each strategy: under the combining strategy, 
all bacteria are under antibiotic stress; under the mixing 
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Figure 2 Emergence of double resistance as a function of the effect of stress on mutation The double resistance emergence for each 

Ms 

strategy is plotted, in log scale, as a function of — , which varies from 1 to 50. Note that lower emergence values indicate higher efficiency of 

the strategies. Mixing is plotted in green, cycling in blue, and combining in red. Double resistance emergence is monotonically increasing in a 
(persistence of antibiotic resistance), so curves with intermediate values of a will be within the boundaries of the curves a = 0 and a = 1 
displayed. The value of a is irrelevant for combining, so only one curve is produced. An intersection of the curves implies that the preference 

Ms 

between two strategies, with respect to the emergence of double resistance, should change. We can see that for a wide range of — , cycling is 

Mr 

the most efficient strategy. When a is low enough, cycling is the most efficient strategy even without SIM. The emergence is taken over 600 
days. The length of each cycle in the cycling strategy is 200 days. Other parameter values: /3 = 0.9, y= 0.03, m = 0.1, X s = 0.1, Xr 1 =0.1, 
Xr = 0.1, t = 0.5. SIM, stressed-induced mutation. 



strategy half of the bacteria of each resistant strain will 
be under antibiotic stress; cycling has the advantage that 
a resistant strain of bacteria will be under antibiotic 
stress only at certain periods of time - when all patients 
are receiving the antibiotic to which this strain is sensi- 
tive. During these periods the hyper-mutating strain 
would typically be at low frequency. This causes the 
increase in double resistance emergence due to SIM to 
be more moderate under cycling than under mixing or 
combining. Persistence of antibiotic resistance, a , also 
has a strong influence on the dynamics: <7 and \i r appear 
in our equations only within a single term a ■ fi r , thus 

Ms 

low values of <7 and high values of — have similar 

effects on the dynamics. The special case where we 
assume no stress induction and a = 0 was previously 
studied by Bonhoeffer et al. [19] and our results are 
consistent with theirs. However, our conclusions hold 
for intermediate values of a as well, assuming mutation 
rates are sufficiently dependent on stress. 

In addition to the parameters pertaining to the stress 
induction mechanisms and within-host selection, the 
parameters determining E SIM also play an important role 



in the emergence of double resistance. To test the 
robustness of our model to changes in these parameters 
we study the emergence of double resistance for 10 
random sets of parameters. The values of m were cho- 
sen from a generalized extreme value distribution fitted 
to the length of patients' stay in observed data [41]. 

Antibiotic clearance rates, —, were chosen from a uni- 

r 

form distribution on [2,14] (deduced from optimal treat- 
ment estimations in [42]). Spontaneous clearance, J, is 
then chosen from a uniform distribution on [0,r] , and 
Infection rates, /3, were chosen from a log uniform dis- 
tribution on [0.001,1] [10]. Four more values were cho- 
sen from the uniform distribution on 0[1], and then 
normalized to determine the entrance frequencies 

^s, / ^R 2 / 1 — ^s — ^.R, — ^R 2 • 

In Figure 3 we plot the value of a moving average (see 
Methods) of the emergence of double resistance, where 
the horizontal axis is taken as the rate of clearance due 
to antibiotic usage (r ). The moving average of double 
resistance under mixing is plotted in green, under 
cycling in blue, and under combining in red. When the 
values of the moving average for mixing and cycling are 
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Figure 3 Robustness of the dynamics describing resistance acquisition by mutation The dynamics describing double resistance 
emergence by mutation were studied for 10 4 random sets of parameters (derived from the distributions described in the main text). The 
moving average of the emergence of double resistance (see methods) is plotted as a function of the rate of clearance due to antibiotic usage 
(r). Note that the same data sets are used for a given value of r under each of the three strategies, generating an association between the 
curves. The moving average of mixing is plotted in green, of cycling in blue, and of combining in red. When the values of the moving average 

Ms 

for mixing and cycling are very close only the blue curve is visible. (A) — = 1, a = 0.1. When o~ is low, cycling decreases the emergence of 



Mr 



Ms 



double resistance better than the other strategies even with no influence of stress on mutation rates. (B) — = 1,0"= 1. When o~ is high and 

jl r 

Ms 

stress has no effect on mutation, combining is more efficient than the other strateqies. (C) — = 10, a = 0.1. At low a, increased influence of 

Mr 

stress on mutation rates results in increased emergence by about an order of magnitude (compare with A), where cycling remains the most 
Ms 

efficient strategy. (D) — = 10, a = 1. At high o~, increase in mutation rates due to stress alters the preferred strategy drastically from combining 
Mr 

to cycling (compare with B). 



very close, only the blue curve is visible. Two values of 
the persistence of antibiotic resistances, are presented: 
a = 0.1(panels A, C) and a = 1 (panels B, D). Cycling 
outperforms the other strategies given even a moderate 
influence of stress on mutation rates, or a low enough 
value of O. As r increases, combining becomes slightly 
more effective (as it reduces single resistance frequency 
efficiently) but the effect is weak. Whenever mutation is 
stress-induced, or a has a low value, combining would 
be the least efficient strategy in terms of inhibiting 
emergence of double resistance. 

Different results are obtained when observing the 
mean proportion of infected patients rather than double 



resistance emergence. Figure 4 shows the mean fre- 
quency of infected patients in the hospital, for the same 
random sets of parameters as in Figure 3 (note that 
under condition CI, any changes in parameters affecting 
stress-induction are irrelevant). We can see that com- 
bining always outperforms cycling and mixing in terms 
of minimizing overall infected patients. When antibiotic 
resistance persistences , is high and mutation is not 
stress-induced, combining reduces both infection and 
emergence more efficiently than the other strategies 
(Figure 3.B), but only by a few percent. However, under 
stress-induced mutation, or low values of (7, there is a 
conflict between minimizing infection and minimizing 



Obolski and Hadany BMC Medicine 2012, 10:89 
http://www.biomedcentral.eom/1 741 -701 5/1 0/89 



Page 7 of 1 3 



0.78r 



.1 076 " 

"5 

CL 

~u 

CD 0.74 
O 
.03 



O 0.72 - 

c 

g 

o 
a. 
o 



0.7r 



C 

cc 

CD 0.68 




).05 



Mixing 
— Cycling 
Combining 



0.1 



0.15 



0.2 



0.25 



0.3 



0.35 



0.4 



0.45 



0.5 



Figure 4 Robustness of the dynamics describing proportion of infected patients The dynamics described in E1 were studied for 10 4 
random sets of parameters (derived from distributions described in the main text). The moving average of the emergence of double resistance 
(see Methods) is plotted as a function of the rate of clearance due to antibiotic usage (r). Note that the same data sets are used for a given 
value of r under each of the three strategies, generating an association between the curves. The moving average of mixing is plotted in green, 
of cycling in blue, and of combining in red. When the values of the moving average for mixing and cycling are very close, only the blue curve is 
visible. We can see that combining is the most efficient strategy in reducing the proportion of infected patients. When the rate of clearance due 
to antibiotics, r, increases, the relative efficiency of combining increases , while the proportion of infected patients decreases. 



emergence of resistance: combining reduces the propor- 
tion of infected patients by about 1% on average, but 
increases emergence of double resistance by more than 
50% (see Figure 3 and Additional file 1 table SI). Simi- 
larly, when a is high and mutation is not stress-induced, 
mixing slightly reduces both infection and emergence in 
comparison with cycling (Figure 3.B and Additional file 
1, table SI). Under stress-induced mutation, or low 
values of <7, there is a conflict: mixing reduces the pro- 
portion of infected patients (by about 0.1% on average), 
but increases emergence. For a large subspace of realis- 
tic parameter sets [43] double resistance is substantially 
decreased under cycling, relative to mixing, under the 
assumption of SIM [See Additional file 1, Figure SI]. 

The clearance rate due to antibiotic usage, r, mildly 
decreases both infection (Figure 4) and resistance emer- 
gence (Figure 3). As the clearance rate increases, bac- 
teria are eliminated more efficiently, thus the window of 
opportunity for infecting bacteria to acquire double 
resistance narrows. 

Stress-induced horizontal gene transfer 

Another mechanism for acquiring antibiotic resistance is 
HGT. Acquiring antibiotic resistance through HGT is a 
process that depends on the rate of encounters between 
bacteria of different resistant strains, and on the prob- 
ability of bacteria to donate and receive genetic material. 
Bacteria of one strain will encounter bacteria of another 



strain at a rate proportional to the amount of interac- 
tions between patients infected by these bacteria. Hence, 
we define 



-J 



r-HGT (U) := R^U^Wdt 



For the rate of encounters between bacteria of differ- 
ent resistant strains we multiply E HGT {U) by a constant 
C, which denotes the rate of bacterial transmission from 
one patient to another (by hospital staff, direct contact, 
and so on). We denote by d the probability that bacteria 
will donate genetic material, and by r the probability 
that bacteria will receive genetic material, given that a 
bacterial transmission event has occurred. The probabil- 
ity of successful HGT between bacteria, given bacterial 
transmission, is thus r ■ d . However, when stress- 
induced HGT is considered, r and d are no longer con- 
stant, as different treatment strategies create varying 
levels of stress for different bacteria. We will define r r , 
d r to be the probabilities of receiving and donating 
genetic material for bacteria which are not under anti- 
biotic stress. Similarly, r s and d s will be the probabilities 
of receiving and donating genetic material for bacteria 
which are under antibiotic stress. For example, if all 
patients are treated with the same drug, bacteria resis- 
tant to that drug would perform HGT with probabilities 
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r r and d r , while bacteria sensitive to it would perform 
HGT with probabilities r s and d s . Stress-induced HGT 
will be expressed by the conditions r s >r r and d s >d r . 

Note that when different patients are treated with dif- 
ferent types of antibiotics, bacteria may be transported 
from a non-stressful environment to a stressful one and 
vice versa, (p will represent the time required for the cel- 
lular mechanisms to induce or repress HGT. Namely, 
when (p = 0 a change in the bacteria's well-being is 
immediately translated to a change in HGT rates, 
whereas for larger values of (p the stress-induction 
mechanisms react more slowly to changes in environ- 
mental stress. We define the extreme case of (p = 1 as 
the case where the HGT rates of bacteria transmitted 
from patient A to patient B depend only on the stress 
the bacteria experienced while residing in patient A. In 
analogy to inequality (CI), we concentrate on the case 
where the influence of non-resistant bacteria on the 
emergence of double resistance is negligible. This can 
be formally expressed as: 



Cd s r s SR «; k R ,m + fiXRi, z = 1,2 



(C2) 



Using the parameters and assumptions described 
above, we denote by Qhgt (U) the emergence of double- 
resistant bacteria through stress-induced HGT under 
treatment strategy U, and derive £,hgt f° r each of the 
strategies [See Additional file 1, section 1, equations 1- 
3]. The combining strategy is the simplest to model, 
since all patients are treated with at least one effective 
antibiotic. Therefore, the emergence of double resistance 
under combining is 

%HGT{combining) = Cr s d s E H c,T (combining). 

Emergence under cycling can be broken into time 
intervals in which only one drug is used. Assuming 
there is an equal amount of such intervals for each drug 
we get 



%HGT{cycling) = C-(r r d s o + r s d r a + r r d s + r4 r )E HGT {q>cling). 

Note that for both combining and for cycling, (p does 
not appear in the term describing double resistance 
emergence, since under both strategies bacteria will not 
experience an environmental change when transported 
from one patient to another. For mixing, the heteroge- 
neity of the environment due to the usage of different 
antibiotics requires the use of (p : 

%HGr{™udng) = C—[r s [d s (p + d,{\ — </>)) + r,d s rj + r s d r + r r (d r <p + d s (l — 0))<x+ 
8 

d s (r s (jt + (1 — 0)r r ) + d s r T + d r r s a + d T {r r (f> + ^(1 — $))(f]EHcr[mixing) 

Analysis of £hgt for each strategy shows that high 

enough values of min {— , — [will cause the ratio 

[d r r r \ 



£hgt [cycling) 

; : to be arbitrarily close to zero [See 

§ hgt (combining) 

Additional file 1, section 1.1, equation (6)]. In other 
words, the more HGT becomes affected by stress in 
both the donor and the recipient ends, the more effi- 
cient cycling becomes relative to combining. This occurs 
due to the fact that when cycling is applied only one 
type of resistant bacteria is stressed at any given time. 
For double resistant bacteria to emerge through HGT, a 
donor and a recipient of different resistance types are 
needed, but under cycling they will never be simulta- 
neously under stress. In contrast, when the combining 
strategy is applied all bacteria are stressed, including 
potential donors and recipients alike. 

The efficiency of the mixing strategy depends on (p. 
This results from the fact that under mixing, different 
patients are treated with different drugs. When (p is low, 
bacteria transported from one patient to another change 
their HGT probabilities according to the drug taken by 
the new host. Therefore, for both mixing and cycling, 
bacterial HGT events which produce double resistance 
will rarely occur when both donor and recipient are 
stressed. In contrast, when (p is high, bacteria trans- 
ported from one patient to another retain HGT prob- 
abilities complying to the antibiotic treatment of their 
former host, allowing the emergence of double resis- 
tance through HGT in which both donor and recipient 
are stressed. 

When exploring the effects of stress induced HGT we 
take — = — , since we have shown that min { — } is 

d T T r [d r r r \ 

the dominant factor in the dynamics [See Additional file 

1, section 1.1 equations 4-6], and in order to avoid 

assumptions about the role of donor and recipient in 

the dynamics. For ease of notation we define 

d s r s . \d s r s 
0 := — = - = mm — , - 
d r r T [ d r ry 

Figure 5 shows the rate of double resistance emer- 
gence for each strategy as a function of 6. We can see 
that the rate of double resistance emergence depends 
also on (p (compare panels A and B of Figure 5), and to 
a lesser extent on c (dashed and solid curves in each 
panel correspond to <7 = 1 and 0 = 0, respectively). 

The ability of a strategy to minimize E HGT plays an 
important role in the emergence of double resistance 
through HGT as well. Since the dynamics determining 
E HGT are complex we will focus only on several impor- 
tant factors. First, we note that when the cycle length 
approaches zero, double-resistance emergence under 
cycling converges to double-resistance emergence under 
mixing, as was shown in [10]. Second, when there is 
strong asymmetry between the frequencies of different 
single resistant bacterial infections, cycling performs 
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Figure 5 Emergence of double resistance as a function of the effect of stress on HGT. Double resistance emergence for each strategy is 
plotted, in log scale, as a function of 8, the effect of stress on HGT. Mixing is plotted in green, cycling in blue, and combining in red. We see 
that stress-induced HGT changes the efficiency of the strategies so that cycling is the most efficient strategy for a wide range of 6 values - note 
that lower emergence indicates a more efficient strategy. Panel A shows results for <p = 1, whereas in panel B <p = 0. We can see that the 
double resistance emergence for mixing indeed increases dramatically for <p = 1 in comparison with (p = 0 , matching the intuitive explanation 
(see text). An intersection of the curves implies that the preference between two strategies should change. The mean is taken over 600 days, 
where the length of each cycle in the cycling strategy is 200 days. Other parameter values: /3 = 0.9, y= 0.03, m = 0.1, X s = 0.1, Xr 1 =0.1, 
Aj} 2 = 0.1, x = 0.5, C = 1 (note that changing C will just rescale the horizontal axis, so its value is arbitrary). HGT, horizontal gene transfer. 



poorly. If the current drug used is ineffective against the 
bacterial strain with the high entrance rates it will allow 
the incidence of that strain to increase rapidly. This was 
shown by Bergstrom et al. [10], who predicted that 
under high asymmetry of single resistance entrance 
rates cycling would not be an efficient strategy. While 
this is true for a model assuming constant HGT, we 
argue that when HGT is stress-induced, and the 
response to environmental changes is not immediate 
(that is, (p > 0 ), cycling can still be the most efficient 
strategy. When the response is immediate, mixing will 
tend to minimize double resistance emergence best. 
Combining will rarely minimize double resistance better 
than the other two strategies when HGT is stress 
induced. These conclusions were obtained by studying 
two prominent factors that influence the emergence of 
double resistance: the ratio of entrance rates [See Addi- 
tional file 1, section 3] and the effect of stress on HGT. 
We have randomly selected 10 4 parameter sets (the 



same parameter sampling as in Figure 3) and for each 
parameter set computed the moving average of the 
emergence of double resistance by HGT for the three 

strategies, as a function of entrance rates' ratio — — (in 

log scale). The moving average of mixing is plotted in 
green, of cycling in blue, and of combining in red. 
When the values of the moving average for mixing and 
cycling are very close, only a blue curve is visible. We 
can see that when HGT is completely unaltered by 
stress, the relative efficiency of the strategies depends on 
the persistence of antibiotic resistances . When O is 
high, combining is the best strategy at minimizing dou- 
ble resistance for most parameter sets (Figure 6.B). 
However, when a decreases, mixing becomes a promi- 
nent strategy as well, with cycling close behind it (Figure 
6. A). In Figure 6.C and 6.D, we allow mild stress 
induced HGT and take 6 = 10. a has a lesser influence 
when HGT is stress-induced [See Additional file 1, 
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Note that the same data sets are used 



emergence of double resistance (see Methods) is plotted as a function of the entrance rates' ratio 

for a given value of under each of the three strategies, generating an association between the curves. The emergence rate under mixing is 

plotted in green, under cycling in blue, and under combining in red. When the emergence rates under mixing and cycling are very close, only 
the blue curve is visible. In (A) and (B), HGT is not affected by stress at all (8 = 1) so the delay of stress-induction mechanisms (<p) is irrelevant. In 
contrast, persistence of antibiotic resistance (a) is influential in such cases, and we present two extreme values. (A) 9 = \p = 0.1. Cycling and 
mixing are the dominant strategies, where their relative efficiency decreases at asymmetric entrance rates. (B) 8 = 1, a = 1 . Combining is the 
preferred strategy for inhibiting double resistance emergence. In (C) and (D) HGT is stress induced with 8 = 10, so while cp is now an influential 
parameter, the effect of a is very minor (Additional file 1, section 1) and we present only the value a = 0.5. (Q 8 = 10, <p = 1, a = 0.5 . Mixing is 
slightly more preferable than cycling, and cycling is most inefficient when entrance rates are asymmetric. (D) 8 = 10, <p = 0, a = 0.5 . Cycling is 
the preferred strategy for inhibiting double resistance emergence. HGT, horizontal gene transfer. 



section 1.1] and, therefore, we set it at an intermediate 
value of <7 = 0.5 . Cycling and mixing now become more 
efficient strategies than combining. We can also see the 
influence of (p on whether cycling or mixing will be the 
most efficient strategy in terms of inhibiting double 
resistance emergence (The influences of other para- 
meters are shown in Figure S2 in Additional file 1). 

As we mentioned above, combining always outper- 
forms cycling and mixing in terms of minimizing overall 
infected patients, but the benefit is of only a few percent 
(Figure 4 and Additional file 1, table SI). When antibiotic 
resistance persistences, is low, even without stress 
induction, emergence of double resistance under com- 
bining is higher by 30% than under other strategies [See 
Additional file 1, table SI]. When stress induction is sig- 
nificant, double resistance emergence under combining is 
more than tenfold higher than under cycling, and more 



than 30% higher than under mixing (Figure 6 and Addi- 
tional file 1, table SI). 

Discussion 

Several conclusions can be derived from our mathemati- 
cal model. We have shown that stress-induced genetic 
variation can have a drastic influence on the emergence 
of double resistance, and should be considered when 
deciding on a hospital wide strategy of antibiotic usage. 
Although always slightly more efficient than other stra- 
tegies in decreasing the incidence of single resistant 
infections, the strategy of combining performs very 
poorly in inhibiting double resistance emergence when 
genetic variation is stress-induced. This holds true 
despite the fact that under the combining strategy all 
patients receive effective treatment, and even though we 
disregard the toxic effects of combining antibiotics for 
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the patient and the economic burden it carries for the 
population [11,44]. 

Cycling is the preferred strategy with respect to the 
acquisition of resistance through SIM. Low persistence of 
antibiotic resistance (a) further amplifies the effects of 
SIM and increases the relative efficiency of cycling. In the 
presence of stress-induced HGT, cycling and mixing are 
the favored strategies, and the preference between them 
is determined by how fast the bacteria respond to envir- 
onmental changes (the parameter (p in our model). If 
changes in HGT frequencies in response to antibiotic 
stress are rapid, mixing is the preferred strategy, whereas 
slow response to stress would tilt the scales in favor of 
cycling. We should note that our predictions hold even 
for a very mild increase of HGT and mutation rates 
under antibiotic stress (Figures 3 and 6) in comparison 
with those described in the literature [25,26,30,31]. 
Higher dependence of variation on stress leads to results 
which are more robust to changes in other parameter 
values. 

There are several criteria which are used to evaluate the 
efficiency of an antibiotic strategy: reduction of total infec- 
tion burden; single resistance minimization; and inhibition 
of multiple resistance emergence [10,19,20,43,45,46]. We 
compared two measures of treatment efficiency: propor- 
tion of infected patients and emergence of double resis- 
tance. The latter is of interest mainly in a population 
where double resistance bacteria are still at a low fre- 
quency, thus we focused on that scenario. There is rarely a 
strategy which is ideal for both infection and emergence of 
double resistance at the same time. A strategy that is suc- 
cessful at reducing infection applies more accurate treat- 
ment, and this has two effects on emergence: on the one 
hand, eliminating infection and thus minimizing the bac- 
teria that would become resistant. On the other hand, 
treatment creates selective pressure and potentially stress 
induced variation - two factors that might lead to a faster 
generation of resistance. We find that combining always 
outperforms mixing and cycling by a small amount, when 
it comes to minimizing infection. When considering 
double resistance emergence, both mixing and cycling out- 
perform combining substantially when variation is stress- 
induced. This contrast should be taken into account when 
deciding on a treatment strategy. 

We make several assumptions that should be discussed 
explicitly. First, we did not consider the possible fitness 
cost of antibiotic resistance. This is consistent with 
recent evidence suggesting that compensatory mechan- 
isms reduce such cost to a low level [27,47,48]. Addition- 
ally, for stress-induced HGT to have substantial influence 
on the dynamics we require that both donor and recipi- 
ent probabilities of HGT would increase with stress. It 
was shown in [30] that acquisition of resistance through 
transformation increases with the recipient's stress. The 



amount of genetic material available for transformation 
in the environment is influenced by the death of bacteria, 
and is therefore dependent on the stress that the donor 
bacteria experience. This is particularly true when phages 
cause lysis of their host when the host is stressed [32,33]. 
Similarly, donation and acquisition of conjugative ele- 
ments were each, separately, shown to increase under 
stress [31,49,50]. Another issue we did not address is the 
influence of stochastic events on the dynamics. The 
population size within a hospital unit might be small 
enough for stochastic events such as epidemic outbursts 
of bacteria and extinction of rare bacterial strains for 
long periods of time, to be very influential [45]. Human 
errors in the form of dosage errors, lack of compliance to 
hospital guidelines and so on can be another source of 
stochastic noise that might shift the dynamics from the 
deterministic expectation described here. 

Furthermore, the values of certain parameters might be 
different for different patients. For instance, elderly 
patients might be more susceptible to bacterial infections 
than other patients [51,52]. This could be expressed by 
modeling compartments of patients with different para- 
meter values (in this example, higher f} values) in accor- 
dance with the epidemiological data. Finally, the relative 
efficiency of the different drugs was assumed to be equal, 
and no drug interactions were considered - two factors 
that may further affect the evolution of resistance [53]. 
Future work could address these matters explicitly. 

Our model points to several directions in which 
empirical data can guide the planning of efficient treat- 
ment strategies. First, it is important to understand 
whether a pathogen acquires resistance primarily through 
mutation or through HGT. Second, it is important to 
estimate the persistence of antibiotic resistant bacteria 
within hosts not currently treated with antibiotics effec- 
tive against those bacteria (the parameter a in our 
model). Finally, we would like to directly assess the 
degree to which stress, and in particular antibiotic stress, 
increases the rates of bacterial mutation and HGT. 
Obtaining such data would be an important step in the 
ongoing struggle against multi-drug resistance. We 
believe that obtaining such precise data will help to 
decrease the prevalence of multiple resistance strains in 
bacterial pathogens which have already shown to increase 
genetic variation under stress, such as P. aeruginosa, S. 
pneumoniae, E. coli and V. cholerae [26,30,31]. 

Conclusions 

In conclusion, our work presents an important factor 
thus far overlooked when planning antibiotic treatment 
strategies, namely the effect of stress on genetic variation. 
We show that considering the effects of stress-induced 
genetic variation alters the results of existing theoretical 
models: specifically, combining antibiotics may result in 
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an increased rate of emergence of double resistant bac- 
teria, whereas cycling antibiotics can be more effective 
than previously thought. Applying our predictions to spe- 
cific pathogens would require better empirical evaluation 
of a few key parameters that affect the dynamics of dou- 
ble resistance emergence. We make specific predictions 
regarding the parameter values that would favor particu- 
lar treatment strategies, suggesting that further investiga- 
tion of stress-induced variation and its mechanisms 
might have crucial importance for combating multiple 
antibiotic resistance. 

Additional material 
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